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In this work we model the quintessence potential in a Taylor series expansion, up to second order, 
around the present-day value of the scalar field. The field is evolved in a thawing regime assuming 
zero initial velocity. We use the latest data from the Planck satellite, baryonic acoustic oscillations 
observations from the Sloan Digital Sky Survey, and Supernovae luminosity distance information 
from Union2.1 to constrain our models parameters, and also include perturbation growth data from 
the WiggleZ, BOSS and the 6dF surveys. The supernova data provide the strongest individual 
constraint on the potential parameters. We show that the growth data performance is competitive 
with the other datasets in constraining the dark energy parameters we introduce. We also conclude 
that the combined constraints we obtain for our model parameters, when compared to previous 
works of nearly a decade ago, have shown only modest improvement, even with new growth of 
structure data added to previously-existent types of data. 

PACS numbers: 98.80.-k arXiv:1501.02678 


I. INTRODUCTION 

Dark energy is one of the great discoveries of the 20th 
century, and remains one of the greatest scientific puz¬ 
zles of the 21st century. While the observational evidence 
does not favour any particular physical mechanism be¬ 
yond a cosmological constant, there are many ideas for 
mechanisms other than the simplest concordance ACDM 
model [T]. The possibility remains that the observed ac¬ 
celeration is caused by a scalar field, quintessence (see 
Refs. [HIS] and references therein). The quintessence field 
affects both the background expansion history of the Uni¬ 
verse, and the rate at which matter over-densities grow. 
A key area in current and future dark energy (DE) ob¬ 
servations is the combination of data on the background 
expansion history and the structure-growth history [1]. 
The relation between these two histories is sensitive to 
the properties of the dark energy mechanism, and could 
allow stronger model constraints than either method sep¬ 
arately. 

In this work, assuming quintessence is a valid de¬ 
scription of observations, we derive new constraints on 
the quintessence self-interaction potential and dynami¬ 
cal evolution, based on the latest observational data on 
baryon acoustic oscillations encoded in the galaxy distri¬ 
bution, the cosmic microwave background, and Type la 
supernovae. We also include galaxy survey data measur¬ 
ing the growth of matter perturbations. Earlier work has 
not included growth-of-structure data (e.g. Refs. [SHE], 
though some recent works have placed constraints on in¬ 
verse power-law quintessence models using new growth- 
of-structure data liiin])- Our work is distinct in that 
it uses a direct modelling of the scalar field potential as 
a Taylor expansion and the field dynamics, rather than 
relying on indirect parameterizations. We restrict our¬ 


selves to the case of a ’thawing’ quintessence m field 
with zero initial velocity at early times. 

The coming decade will see a series of observational 
projects, such as the Dark Energy Survey [T^, BOSS 
[I3], BigBOSS [H], WFIRST [15], Planck [l6], HET- 
DEX [T7], DESI [TH], Euclid [121 and LSST [20], all 
of which will probe the growth of structure by means 
of different techniques. Some of the prospects for con¬ 
straining quintessence models with these are discussed in 
Refs. [UlllIlT]. 

II. FORMALISM 
A. Cosmological model 

We follow the approach pioneered in Refs. 016]. We 
briefly review the set-up here, which is fairly straight¬ 
forward. We assume that the quintessence field </> has a 
self-interaction potential V{(f ), that we expand in a se¬ 
ries about the present value of the field which is taken 
(without loss of generality) to be zero. The quintessence 
field obeys the equation 

with the Hubble parameter H given by the Friedmann 
equation 

H = {p + P 4 ,) ■ (2) 

Here p represents all the usual material components: 
dark matter, baryons, neutrinos, and radiation. We will 
use pra to give the matter density (dark matter plus 
baryons) and = 0^/2-|-U(0) the quintessence density. 
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We assume spatial flatness throughout (as motivated by 
cosmic microwave background (CMB) measurements and 
the inflationary paradigm), though the generalization to 
the non-flat case would be straightforward. Our compu¬ 
tations start at a redshift of one million. 

An important quantity, which determines the cosmo¬ 
logical effects we consider from the quintessence field, is 
the equation of state 


_ <^V2 - ¥{(!,) 

P4> <^2/2 + V"(0) ■ 


In order to have close to —1, mimicking ACDM, 
the evolution of the held should be potential dominated, 
which is usually called the slow-roll regime. We focus on 
thawing quintessence, where the field is initially frozen 
due to high damping from the Hubble friction term, 
meaning its energy density is fairly constant at early 
times. However, at late times, the held starts thawing 
and slowly rolling on its potential towards larger values 
of [2] . 


B. Parameterizations and priors 

1. Power series 

We will use a Taylor expansion power series to model 
the potential as 

Vi^) = Vo + Vi(l) + V2^'^ + ..., (4) 

and work in units of the reduced Planck Mass, Mp = 
I/VSttG, defining (j) in dimensionless units. To have vi¬ 
able flat models, we set Vq = P^/Pcri = which will 
be derived from the present-day energy density of mat¬ 
ter, Hni, due to the spatial flatness condition. Note that 
is the sum of the present-day cold dark matter den¬ 
sity and baryon density divided by the present-day crit¬ 
ical density, Qc and Hb respectively. also includes 
the contribution of massive neutrinos, for which we fix 
= 0.06 eV. 

For all our models we have set 0 = 0 at the initial 
redshift, Zi « 10®, so that the scalar held starts at rest. 
We perform a simple binary search to determine the right 
= 0i that allows us to recover a flat cosmology 
today. We then rescale the Taylor series parameters of 
the potential by the value of (j){z = 0) = (()f so that the 
present-day value of (j) is set to zero. As an example, for 
the first-order series expansion of V ((/>), this means that 

Ho =Vo + Vic/>f, 

lA =Hi, 

(^i = (jji- (/){ 

and so on for higher-order expansions of the potential. 
Hence, Vq will correspond to the present-day value of the 
quintessence potential. Effectively, we will be constrain¬ 
ing the rescaled parameters, which we indicate with the 
tilde. 


2. Priors 

In this work, as detailed before, we expand the 
quintessence potential in a Taylor series expansion. We 
limit ourselves to a second-order expansion, which means 
we will have three dark energy parameters describing our 
model, apart from the cosmological parameters such as 
ujc and ln(lO^°As). For these parameters we have kept 
the same prior ranges used in the Planck analysis [22) . 
which are default for CosmoMC [53]. We treat CMB 
foreground parameters identically to Ref. [55] . 

In the previous section we defined the zeroth-order 
term as Vq. This is constrained by to guarantee 
models with a flat cosmology and when only this term 
is present it corresponds to the present-day energy den¬ 
sity of our dark energy component. If the potential is not 
flat there will be an additional contribution to the present 
field energy density from its kinetic energy, which data 
will constrain to be fairly small corresponding to close 
to —1. 

We have imposed a flat prior on Vi in the range [—4,4] 
for the first order expansion of the potential. The reason 
for this is that this range includes, conservatively, the 
one considered in previous similar works [3IS]. Given 
the time elapsed since those works were done, we expect 
that our range won’t have a direct impact on the results 
and should be quite constrained by the data. 

Lastly, there is the second-order term, The previ¬ 
ous works we refer to were not able to detect significant 
constraints on this parameter. After some inspection, we 
believe our results for this parameter are also very prior 
related, for reasons we will discuss later. Hence, without 
any particular motivation, we have set a flat prior ranging 
between [—10,10]. For the Vi parameter, we reduced the 
prior to [—2, 2], in order to avoid fast rolling of the scalar 
field for large values of V 2 , which could jeopardize the 
success of the numerical evaluation of the cosmological 
evolution predicted by our models. 


III. DATA ANALYSIS 
A. Observables 

The analysis done in this work was performed using 
the publicly-available CosmoMC code [53]. For that, we 
modified the corresponding quintessence module. For our 
analysis, we have considered the first Planck data release 
[55] plus the WMAP polarization data [53]. We have in¬ 
cluded CMB lensing in our analysis, as well as the SPT 
[53] and ACT [53] data. We have also included the bary- 
onic acoustic oscillations data from the seven-year and 
nine-year release data sets from the Sloan Digital Sky 
Survey (SDSS) [57] [5B]. Furthermore, we use the Union 
2.1 580 SNIa catalogue from the Supernova Cosmology 
Project (SCP) [29], where we conservatively include sys¬ 
tematic errors and marginalize over TLg as detailed in the 
appendix of Ref. [30] . 
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Lastly, we include growth of structure data from the 
WiggleZ Dark Energy Survey [3T], the BOSS survey (32] 
and the 6dF Galaxy Survey [33] • The two-dimensional 
power spectrum data from these surveys were used to fit 
for both the redshift-space distortion effect (which mea¬ 
sures the rate of growth of structure /erg) and the Alcock- 
Paezynski distortion of the survey geometry. The mod¬ 
ule we use compares our predictions against both types 
of measurement. Further details of these measurements 
can be found in Ref. [34] . 

B. Parameter estimation 

The parameter space we study will be 

© = {l006iMC,T,Wc,Wb,ln(l0^°As) ,ns, Vi,V2} , (5) 

where the cosmological parameters have the same mean¬ 
ing as in Ref. [33|. The parameter estimation is carried 
out using an MCMC approach. The posterior probability 
of the parameters 0, given the data and a prior proba¬ 
bility distribution 11(0), is 

P(©|data) = 5^1e-(xsNia+xcMB+XBAo+XGRo)/2 ^ ( 5 ) 

where the is a measure of the goodness of 
fit between the model theoretical predictions and 
the observed values of a certain physical quantity; 
Z = J £(data|©)n(©)d© is a normalization constant, 
irrelevant for parameter fitting. 

IV. RESULTS 
A. Parameter estimation 

1. Cosmological constant 

In this subsection, we show the marginalized 2-d con¬ 
tour of the Planck constraints assuming a pure cosmo¬ 
logical constant in our model. Therefore, we have set 
hi = 0 for i ^ 0. Hence, Vb = Vo is just a derived pa¬ 
rameter defined as the present-day energy density of the 
scalar field, with a value set to yield a flat cosmology. In 
Fig.jT] we can see that our results match those from the 
Planck data release [22], as expected. 

2. Linear potential 

Here, we present the results for the linearized 
quintessence potential model, for which we have set 
Vb = 0. In Fig. [^ we show the marginalized 2-d con¬ 
tours as well as the one-dimensional marginalized prob¬ 
ability distributions for all the parameters considered in 
this model. We show the constraints from the individual 
data sets as well from their combination. 



FIG. 1: The 2-d marginalized Planek constraints over the 
zeroth-order Taylor series potential model, representing an 
effective cosmological constant. We recover the results ob¬ 
tained by the Planek analysis for a flat Universe. 

In the top plot of Fig. [^ we have the constraints placed 
by the different data sets, individually, on the cosmologi¬ 
cal parameters, PIq and Hm, and the dark energy parame¬ 
ters, Vq and Ui = Vi. The constraints on Vq closely follow 
those obtained for Dm since, as specified before, Vq is ac¬ 
tually a derived parameter, constrained from imposing 
a flat cosmology in our models. The larger uncertainty 
in Vb is due to the inclusion of an additional parameter 
introducing new degeneracies in explaining the data. We 
note that the Planck data are clearly superior in con¬ 
straining the traditional cosmological parameters (here 
PIq and Hm), but that this superiority does not extend 
to the dark energy parameters. 

The supernova data are the most effective in con¬ 
straining the DE parameters, with the smallest limits on 
the posterior probability distribution of V/. The other 
datasets give comparable performance to each other. 

Combining the four data sets, as seen on the lower 
plot of Fig. [^ the confidence limits improve consider¬ 
ably. The confidence limits resulting from the four data 
sets combined can be seen in Table [l] Nevertheless, the 
improvement is modest compared to previous works of 
Refs. [3(3], with most of the constraining power coming 
from the Supernovae alone and the contributions from 
the other data sets are not particularly enhancing. As 
expected the posterior probability distribution of Vi is 
symmetric around zero, the perfect cosmological constant 
case, which is well within the 68% confidence limits. 

It is interesting to see how the individual data sets 
behave and how they combine, particularly in the 
plane. We observe that the constraints from Supernovae 
and Planck exhibit opposite trends in how the preferred 
|Vi| changes with In the Supernovae case, a tilted 
potential compensates a smaller matter energy density in 
the luminosity distance because a higher tilt corresponds. 
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TABLE I: 68% confidence limits of the first and second-order potential model parameters from the different data sets considered 
and their combination. We present the 95% upper limit on the absolute value of Vi, given its reflection symmetry. We don’t 
show confidence ranges for V 2 as the models at the prior edges are not considerably disfavored compared to the best-fit. 




Ho 

Om 

Vb 

|ld| 

V 2 (best-fit) 


Planck 

61 ±6 

0.39 ±0.09 

0.4 ±0.3 

|Ui| < 3.3 



Union 

78± 15 

0.26 ±0.06 

0.67 ±0.05 

|Ui| < 1.3 


Linear 

BAO 

79± 13 

0.29 ±0.12 

0.5 ±0.3 

|Ui| <3.2 



Growth 

85± 11 

0.26 ±0.06 

0.41 ±0.25 

|Ui| < 3.3 



Combined 

67.2 ± 1.0 

0.311 ±0.011 

0.66 ±0.03 

|i4| < 0.77 



Planck 

66 ±3 

0.33 ±0.04 

0.60 ±0.10 

|Ui| < 2.8 



Union 

79± 15 

0.27 ±0.05 

0.68 ±0.06 

|Ui| < 2.5 


Quadratic 

BAO 

82± 13 

0.25 ±0.10 

0.68 ±0.13 

|Ui| <2.7 



Growth 

80± 12 

0.27 ±0.06 

0.64 ±0.11 

|Ui| < 2.8 



Combined 

67.2 ± 1.1 

0.311 ±0.012 

0.64 ±0.06 

|Ui| <2.7 
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effectively, to a higher dark energy equation of state (i.e. 
higher than —1). Hence, the dark energy contribution at 
higher redshift increases, compensating for a smaller Hm, 
leading to the left-sided boomerang shape we observe. 

As for the Planck constraints, the distance-related ob¬ 
servables considered are the acoustic scale, I a, and the 
shift parameter R [35]. I a corresponds to the ratio be¬ 
tween the angular-diameter distance. Da, and the co¬ 
moving sound horizon, rg, at photon decoupling. Increas¬ 
ing decreases both quantities, and a tilted potential 
also decreases Da due to a higher dark energy effective 
equation of state, compensating for higher matter energy 
densities. Lastly, R is the product between Hm and Da- 
Increasing the matter density is, therefore, compensated 
by a higher Vi. All these effects end up contributing for 
the right-shaped boomerang contour we observe. 

A further aspect worth discussing is the clear bimodal¬ 
ity in the Hq-Vi contours produced by the growth data 
(also present in the Hm-Hi contours) in Fig. As an 
increase in Hq increases erg by altering the shape of the 
power spectrum, a higher density of dark energy is re¬ 
quired in order to remain consistent with the growth 
measurements of /ergj through lowering /. 

Finally, in Fig. we see the constraints obtained in 
the WQ-Wa plane, where dw/da = —Wa was computed 
at z = 0. We note how our results restrict themselves 
to a thawing regime, with the scalar field becoming free 
at later times to roll down the potential and increasing 
the effective equation of state when Vi 7 ^ 0. Hence our 
constraints are restricted to the icq > and Wa < 0 , 
with the data preferring the ACDM case, as expected. 

3. Quadratic potential 

In this subsection, we show the results for the quadratic 
quintessence potential model, where we vary both Vi and 
V 2 . In Fig. 1^ we have the marginalized 2-d contours as 
well as the one-dimensional marginalized probability dis¬ 


tributions of all the parameters considered. We show the 
constraints from the individual data sets as well from 
their combination. The confidence limits of the parame¬ 
ters that define this model can be seen in Table [T| 

Inevitably the constraints over our parameters get 
worse by considering a higher order of the potential. As 
in the previous case, the Supernovae data performs bet¬ 
ter (although now only marginally better) in constraining 
the dark energy parameters over the other three available 
data sets. The confidence limits in Vo and Vi almost dou¬ 
bled compared to the first-order potential counterparts. 
The limits on the 1-d posterior probability distribution of 
V 2 extend to the prior we have defined for this parameter 
on all the data sets considered, in a clear indication of the 
difficulty of the data in constraining it. Our understand¬ 
ing is that the limits would always be cut by the prior we 
would define, had it been smaller or larger. For instance, 
the positive, upper limit would extend indefinitely, as the 
data will always accept models with a large V 2 as long 
as Vi remains zero. This is because, in our search for 
the initial conditions, that means the field will remain 
immobile at the minimum (cj) = 0 ) of the convex poten¬ 
tial, producing, effectively, a ACDM-like model, which is 
in strong agreement with the data. Hence, one obtains 
the narrow region extending to high values of V 2 in the 
2 -dimensional Vi — V 2 contours. 

As for the bifurcated shape in the Vi — V 2 contours 
when considering the negative region of V 2 , particularly 
for smaller, negative values of it, is due to the fact that 
it is increasingly difficult to obtain a zero-valued Vi in 
these conditions. The concave nature of the potential 
makes it hard for the field to remain still ai (f) = 0, as the 
smallest tilt will make the field roll considerably. Hence, 
the displacement of the field will always produce a non¬ 
zero value of Vi, creating the shape we can see in Fig. 
This could extend indefinitely according to the negative 
prior limit on V 2 , with the increasing displacement of the 
field producing a widening of the legged shape figure. 

In Fig. we have the constraining contours on the 
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FIG. 2: The 2-d contours and one-dimensional probability distributions for the first-order Taylor series potential model 
parameters. The top plot shows the constraints from the individual data sets (68% and 95% contours), while the bottom figure 
presents the constraints when all of the four data sets are combined (68% and 95% contours). The Planck contours include 
Planck, ACT and SPT data, as well as CMB leasing and WMAP polarization. 
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iDq dw/da Vi 


FIG. 3: The 2-d contours in the Wo-Wa plane and the re¬ 
spective marginalized posterior probabilities distributions for 
the first (upper plot) and second-order (bottom plot) Taylor 
series potential, dw/da = —Wa was computed at z = 0. Since 
the x-axis limits are so different between the different orders 
of the potential expansion, we do not join both plots. 


Wq - Wa plane for the second-order Taylor expanded po¬ 
tential. We observe something very similar to the lin¬ 
ear potential results, with our results constrained to a 
thawing regime. We do have now a very small region 
of Wa > 0, which is however not indicative of a tradi¬ 
tional freezing regime. This actually is a result of the 


field rolling up the potential close to z = 0 after passing 
its minimum and being slowed down by the slope of the 
potential. This type of behavior is observed, for instance, 
in the models considered in Ref. [35] ■ We also note the 
very different WQ-Wa ranges between the two models, as 
the quadratic potential is able to cover much more of this 
plane with observationally-permitted models. 

Lastly, we also observe an odd feature in the V 2 -Hq 
contour plots of Fig. For V 2 > 0, we note that the 
allowed values of Vi get significantly smaller. Hence, this 
shrinks Vq and towards larger and smaller values re¬ 
spectively, with affecting the allowed values of Hq 
(towards larger values) due to the known degeneracy be¬ 
tween these two parameters. The same reasoning can be 
applied for V 2 ^ 0: Vi can assume larger values, widening 
Vq and towards smaller and larger values respectively, 
making Hq tend to smaller values and producing a bump 
in the V 2 -Hq contour plot. 


V. CONCLUSIONS 

In this work, we expanded the quintessence scalar field 
potential in a Taylor series around the field’s present-day 
value, which we set to zero by a rescaling of the potential 
expansion parameters. We also limit ourselves to a thaw¬ 
ing scenario by setting the field’s initial velocity to zero. 
We then used CosmoMC to constrain these parameters 
with the latest available data. 

The main addition in this study was the inclusion of 
growth of structure data, which constrains the models 
against observations of redshift distortions of objects and 
measurements of the normalized growth rate, fa^{z). We 
conclude that the growth data performance in constrain¬ 
ing the dark energy parameters we have considered is 
quite competitive compared to the other three data sets 
we have used, with the supernova data providing the 
tightest constraints on the potential parameters. 

We have also shown explicitly how the different data 
sets constrain our model parameters. It is particularly 
interesting how the Supernovae and the Planck contours 
behave in opposite manner in the — Vi plane. This 
could be promising for future surveys, proving to be the 
ideal combination for constraining dynamical dark en¬ 
ergy models due to the observables involved. 

The improvement in constraints we observe is mod¬ 
est relative to our previous works EE], where the data 
used was the WMAP 3-year distance information data 
the SNLS type la Supernova dataset [35], and the 
SDSS and 2dF galaxy redshift survey baryon acoustic 
peak information isaiioi- It is true that the model we 
consider is different by construction, since we restrict our¬ 
selves to thawing quintessence models, but nonetheless 
the constraints on the parameters of the expanded po¬ 
tential did not improve significantly in the many years 
since those works despite the accumulated data of the 
past decade, both of existing types and new growth of 
structure data. However it is notable that those early 
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FIG. 4: As Fig. for the second-order Taylor series potential model parameters. 













































works had a fairly narrow prior on Vi which contributed 
to the quoted strength of the constraint, whereas our 
wider prior shows this parameter to now be well con¬ 
strained by data in the linear potential case. 

A comparable example is the marginal improvement 
of constraints on the DE equation of state, w, or on the 
2 -dimensional contours on the WQ-Wa plane (we suggest 
comparing the results from the WMAP 9-year data with 
the Planck results [55J [23] ). It has been shown that the 
WQ-Wa parameterization is a good representation of thaw¬ 
ing quintessence models |41] . Therefore, given the expec¬ 
tation that next-generation experiments will be unable to 
distinguish thawing quintessence with wq < —0.9 from 
ACDM at the 95% confidence level [42], we predict that 
improvement on the constraints we show in this work to 
be challenging also. 

Despite this, we conclude that for the first-order Taylor 
series expansion of the potential, ACDM is well within 
the preference of the data, with the confidence limits on 
the parameters we consider being Vq = 0.66 ± 0.03 and 
|Ei| < 0.77 (95%). These constraints get significantly 
worse when we consider an extra order on the expansion 
of the potential. The second-order term, V 2 , is uncon¬ 


strained by the data, with the posterior limits on this 
parameter extending to the prior range we define for it, 
possibly indefinitely. 
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